Sequence-dependent spin-selective tunneling along double-stranded DNA 
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We report spin-selective tunneling of electrons along natural and artificial double-stranded DNA 
(dsDNA) sandwiched by nonmagnetic leads. The results reveal that the spin polarization strongly 
depends on the dsDNA sequence and is dominated by its end segment. Both genomic and artificial 
dsDNA could be efficient spin filters. The spin-filtering effects are sensitive to point mutation which 
occurs in the end segment. These results are in good agreement with recent experiments and are 
robust against various types of disorder, and could help for designing DNA-based spintronic devices. 
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The charge transport along DNA molecule has received 
significant attention from scientific researchers over the 
past two decades. [ll-Q In addition to electric charges, the 
DNA molecule could be also used to manipulate the elec- 
tron spin. It was reported that self-assembled monolay- 
ers of double-stranded DNA (dsDNA) can discriminate 
the spin of photoelectrons.0, H| These electrons trans- 
mitted through the dsDNA monolayers are highly polar- 
ized at room temperature and the spin-filtering effects 
are enhanced with increasing the DNA length. |5fl More- 
over, it was demonstrated that even single dsDNA could 
be efficient spin filter. Q The underlying physical mech- 
anism arises from the combination of the dephasing, the 
SO coupling, and the chirality of the DNA molecule. [7[ 
However, the spin polarization vanishes if the dsDNA 
was changed into single-stranded DNA or damaged by 
ultraviolet light. [1-01 

The nitrogenous bases guanine (G), adenine (A), cy- 
tosine (C), and thymine (T), which are four basic ingre- 
dients of the DNA molecule, can constitute thousands of 
various sequences. While natural DNA molecule can be 
extracted from the cells of all living organisms, the arti- 
ficial one could be synthesized in any desired sequence. 
It was shown that the DNA molecule with different se- 
quences could present any transport behavior of con- 



ducting, semiconducting, and insulating, [ctillj One may 
thus expect that different dsDNA would display diverse 
spin-filtering effects. Indeed, the study of spin transport 
along various dsDNA will provide valuable information to 
the physical mechanism and the biological processes, and 
opens up its potential applications in molecular spintron- 
ics. In this Letter, we explore spin-selective tunneling of 
electrons through the dsDNA connected by normal-metal 
leads. Based on a model Hamiltonian which includes the 
SO coupling and the dephasing, the conductance and the 
spin polarization are calculated for a variety of dsDNA. 
Here, the DNA molecules involve genomic and artificial 
ones as well as those employed in the experiments. 0, Q 
The sequences of several typical DNA samples are listed 
in Table HI The genomic dsDNA is extracted from the 
sequence of human chromosome 22 (chr22),[l2j while the 
artificial dsDNA is taken as random sequence and sub- 



stitutional one, e.g., Nickel mean (nm), Copper mean 
(cm), and Triadic Cantor (tc).[l3| All of the substitu- 
tional DNA sequences are constructed by initiating from 
one seed and following a substitution rule. For instance, 
the nml sequence is formed by adopting base A as the 
seed and the substitution rule A— sAGGG, G— >A. 

From the study of numerous dsDNA, we find that the 
spin filtration efficiency presents strong dependence on 
the DNA sequence and is mainly determined by the end 
segment with several base-pairs. Both chr22-based and 
random dsDNA could be very efficient spin filters, while 
the substitutional one exhibits large spin polarization and 
conductance. Besides, the spin-filtering effects arc sen- 
sitive to point mutation which takes place in the end 
segment of the dsDNA. The high spin polarization still 
holds even under the environment-induced on-site energy 
disorder and twist angle disorder. These results could be 
beneficial for building up DNA-based spintronic devices. 

The spin transport along the dsDNA can be described 
by the Hamiltonian: V. = Hbna + H so + Hd + "Hioad + 
%c@, Q where "Hdna = X^=i(X^=i £ jnc] n c jn + 

J2n=i t jnc] n Cjn+i) + J2n=i A„c|„c 2 „ + H.c. is the Hamil- 
tonian of two-leg ladder model, with n the base-pair in- 
dex, j the strand label, and N the DNA length. cJ n = 

fct. , cj„ I ) is the creation operator, e ln is the on-site en- 
ergy, tj n is the intrachain hopping integral, and A„ is the 
interchain hybridization interaction. The second term 

?4o = Sj n {^soc]ri cr n )c j'n+i + H.c.} is the SO Hamilto- 
nian, which stems from the double helix distribution of 
the electrostatic potential of the dsDNA.Q t so is the SO 
coupling strength and 0$ — [a x (sin (pj n + sin<y9 jn+ i) — 
a y (cosip jn + cos^jn+i)] sin^n + 2a z cosO jn , with a x . y . z 
the Pauli matrices, ipj n the cylindrical coordinate of the 
base, and 9j n the helix angle between base n and n + 1 
in the jth strand. In equilibrium position of the dsDNA, 
Wjn = { n ~ and 0j n = with Atp the twist angle. 

The third one U d = T^jnki E ^kb] nk b jnk +t d b\ nk c jn +R.c.) 
is the Hamiltonian of the Buttiker's virtual leads and 
their coupling with each base of the dsDNA, simulat- 
ing the phase-breaking processes due to the inelastic 
scattering with phonons and counterions. 15, Gil The 
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TABLE I: The sequences of the DNA molecules. Here, only the sequence along one strand is presented, while the other can 
be derived according to Watson-Crick base-pairing rules: G pairs with C, and A with T. The first three terms are the DNA 
molecules adopted in the experiments, rdl, rd2, and rd3 are the random sequences, hcl, hc2, and hc3 are the chr22-based 
sequences, and the last four terms are the substitutional ones. 
Name DNA sequence 

sq-26 TTTGTTTGTTTGTTTGTTTTTTTTTT 

sq-40 TCTCAAGAATCGGCATTAGCTCAACTGTCAACTCCTCTTT 

sq-50 TACTCTACCTTCTCAAGAATCGGCATTAGCTCAACTGTCAACTCCTCTTT 

rdl CAATGCAGTCTATCCACCTGACGGACCCCGACCCGCCTTT 

rd2 CAATGCAGTCTATCCACCTGACGGACCCCGACCCGGCTTT 

rd3 CAATGCAGTCTATCCACCTGACGGACCCCGACCCGCCATT 

hcl TAAATAAATAAATAAATAAATAAAATAAATAAAAGCCTTT 

hc2 GGGCCCTGAGGCATGGGCCCAGAAGCATTCCTGTCCCCTT 

hc3 AGCTGGGGAGCAGGGCTCCACTCTGGGAGGGGGGCAGCCT 

nml AGGGAAAAGGGAGGGAGGGAGGGAAAAGGGAAAAGGGAAA 

nm2 ATTTAAAATTTATTTATTTATTTAAAATTTAAAATTTAAA 

cm GAAGGGAAGAAGAAGGGAAGGGAAGGGAAGAAGAAGGGAA 

tc GAGAAAGAGAAAAAAAAAGAGAAAGAGAAAAAAAAAAAAA 



last two terms %i cad + H c = J2k,t3=LM £ Pkal k a i3k + 

^jk(tLa\} ; Cji+t]ia} Rk CjN+^i-C-) represent the real leads, 
and the coupling between these leads and the dsDNA, 
respectively. Finally, the conductances for spin-up (Gf) 
and spin-down (G±) electrons can be calculated by using 
the Landauer-Buttiker formula. Q The spin polarization 
is P s = (Gf — Gi)/(Gf+Gi). Since the current is flowing 
from the left real lead to the right one, the terminal of the 
dsDNA attached to the former is named the beginning, 
while the other terminal is called the end. 

For the dsDNA, Ej n is chosen as the ionization po- 
tential with eq = 8.3, ea = 8.5, Ec — 8.9, and 
ct = 9.0, tj n between identical neighboring bases is 
taken as £gg = 0.11, <aa = 0.22, tec — —0.05, and 
t TT = —0.14, and A„ = —0.3. These parameters are 
extracted from the experimental results 17 - 2(J and first- 
principles calculations 2114271] with the unit eV. tj n be- 
tween different neighboring bases X and Y is set to 
ixY = (£xx + £yy)/2, in accordance with first-principles 
results. 23|-26| The helix angle and the twist one are 
6 = 0.66 rad and Aip = |. The SO coupling is esti- 
mated to t so = 0.01. For the real leads, the parameters 
= = 1 are fixed, while for the virtual ones, the 
dephasing strength is Td — 0.006. 

It was reported that the ionization potential of the 
base is affected sig nificantly by both counterions[2^, [29[ 
and hydration. [30(-|32j Consequently, the environmental 
effects can be properly considered by varying the on-site 
energies. A random variable Wj n is added in each £j n 
to simulate the stochastic population of these counteri- 
ons and water molecules around the dsDNA, with Wj n 
uniformly distributed within the range [—■§-, ■§-] and W 
the disorder degree. Fig. QJa) shows the spin polarization 
P s of poly(A)-poly(T) under the on-site energy disorder, 
as a function of the energy E. It clearly appears that 
P s is large for homogeneous poly(A)-poly(T) and is suf- 
ficiently robust against the on-site energy disorder. This 
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FIG. 1: Energy-dependent P 3 for poly(A)-poly(T) under the 
on-site energy disorder with degree W (a) and of the twist 
angle disorder with degree D (b). (Gf) and (Ps) vs W (c) 
and vs D (d). The inset of (c) shows (Gf) in a wider range of 
W and the dependence can be fitted by the function (G^) oc 
l0~ aW with a = 5.80 ± 0.09 (cyan line). (G t ) and (P s ) are 
averaged in the energy region [9.04,9.32]. All of the results 
are performed for single disorder configuration with N = 40. 
Here, Go = e 2 /h is the quantum conductance. 



can be further demonstrated in Fig. [He), where the av- 
eraged spin polarization (P s ) is shown. One notices that 
(P s ) fluctuates around its equilibrium value of 5.0% at 
W = and the oscillation amplitude is enhanced by W . 
Furthermore, a new energy region of high P s becomes 
more distinct in the case of larger W [see the curves of 
W = 0.16 and 0.3 in Fig.^a)]. On the other hand, the 
averaged conductance (Gf) is decreased by increasing W 
as expected, due to the disorder-induced Anderson local- 
ization effect. The curve of (Gf)-W can be fitted well by 
a simple function (Gf) cx \0~ aW [see inset of Fig. [TJc)]. 

Besides the on-site energy disorder, each base will wa- 
ver around its equilibrium position at finite temperature. 
In this situation, it is reasonable to plus a random vari- 



able dj n in each >pj n , with dj n distributed in the region 
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FIG. 2: Energy-dependence of G^-, Gj., and P s for the ds- 
DNA used in the experiments, (a) G-\ and Gj. for the sq-26 
sequence, (b) Gf for the sq-40 and sq-50 sequences, (c) P a for 
all three dsDNA. The inset of (c) displays P s with F d = 0.01. 



[—|,y] and D the disorder degree. By considering con- 
stant radius R of the dsDNA and arc length l a between 
successive bases to account for the rigid sugar-phosphate 
backbone, @, [33| the helix angle 6j n will be modulated 
according to l a cos9j n = Rifjn+i — <pjn) and the fluc- 
tuations are disregarded in the intrachain hopping inte- 
gral as a first approximation. [2F 



34 |35| It can be seen 



from Fig. QJb) that the curves of P s -E are superposed 
with each other in the context of the twist angle disorder 
only. Accordingly, no fluctuations could be observed in 
the curve of (P s )-D [Fig. QJd)]. Besides, (G t ) will not 
be changed with D, because the SO coupling is much 
smaller than the hopping integral. Therefore, poly(A)- 
poly(T) remains an efficient spin filter even under the 
on-site energy disorder and the twist angle disorder. 

Then we investigate the spin transport through ape- 
riodic dsDNA in the absence of external environment- 
induced disorder. Our results still hold if this disorder is 
included. Let us first discuss the spin transport proper- 
ties of the dsDNA used in the experiments. @,@] Fig.[^a) 
displays the conductances of spin-up (G|) and spin-down 
(Gj.) electrons for sq-26 sequence, while Fig. [2f b) plots 
G^ for sq-40 and sq-50 sequences. As compared with 
homogeneous dsDNA,0 the energy spectrum of aperi- 
odic dsDNA is also separated into the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied 
molecular orbital (LUMO). The conductance is declined 
by increasing the DNA length, because the electrons ex- 
perience stronger scattering in longer dsDNA. 

In addition, one can see from Fig. EJa) that the dis- 
crepancy between G-f- and G^ is more distinct for the 
LUMO band than the HOMO one. Thus, P s is larger 
in the former band than the latter one [Fig. [He)]. More- 
over, P s at E = 9.1 f is, respectively, 9.6%, 38%, and 38% 
for the sq-26, sq-40, and sq-50 sequences, in good agree- 
ment with the experiment. [5j In fact, the obtained P s is 
also consistent with the experimental results by adopting 
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FIG. 3: Energy-dependent P„ for the random dsDNA (a) and 
for the chr22-based one (b). (c) Distribution of P s at E — 9.11 
for various random dsDNA with large P a (right part) and with 
small P 3 (left part) as a comparison. The results are extracted 
from 10 DNA samples. Here, only the end segment in the 
first strand is shown (two sides) and can be obtained for the 
second one according to the base-pairing rules, (d) (G^) and 
(P s ) vs mutation position L for the rdl sequence. 



different from the region [0.003,0.01], e.g., see P s of 
Td = 0.01 in the inset of Fig. [UJc). Besides, although the 
conductances between the sq-40 and sq-50 sequences are 
very different, their spin polarizations are almost identi- 
cal and the difference between the two P s is within 10~ 6 
range, due to the fact that the sq-40 sequence is the end 
part of the sq-50 sequence. This suggests that the spin 
filtration efficiency of the dsDNA is mainly controlled by 
its end segment, which will be substantiated below. 

Next we turn to study the spin polarization of the ran- 
dom and chr22-based dsDNA. Figs, [^a) and Mb) plot 
P s vs E for several typical random and chr22-based se- 
quences, respectively. It is clear from the curves of rdl 
and hci that both random and chr22-based sequences 
could be very efficient spin filters with P s achieving 40%. 
From a statistical study of numerous dsDNA with ex- 
tremely high P s , it reveals that these sequences are termi- 
nated by the segment "CCTTT / GGAAA" in their ends 
[Fig. [UJc)]. We emphasize that all of the investigated ds- 
DNA with TV = 40 will exhibit very high P s around 40% if 
their end segments are replaced by "CCTTT/GGAAA" . 
Besides, the dsDNA could be also very efficient spin filter 
if it has other end segments, as shown in Fig. EHc), where 
a distribution of P s at E = 9.11 is displayed for different 
random dsDNA with 16 end segments. It clearly appears 
that P s is always large for these dsDNA [see the right part 
in Fig. [He)], although P s will vary in a finite range. The 
dsDNA remains efficient spin filter if it is ended by the 
moiety "(C) m TT/(G) m AA" with m the integer (see the 
curve of hc2). However, P s can be dramatically reduced 
by altering the end segment, even if its last base- pair is 
changed [see the left part in Fig. (SJc)]. These are due to 
the fact that the charge will gradually lose its phase and 
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FIG. 4: Distribution functions of P s for the random and 

chr22-based dsDNA at E = 9.11. The inset shows the corre- _ T _ _„ . , _ _ , ,-, r ■ i ,• 

sponding statistics of P s at E = 9.03. Here, N = 40. FIG ' 5: , Energy-dependence of G t and P s for several substi- 

tutional sequences of DNA molecules. 



spin memory while transmitting along the dsDNA. The 
longer distance the charge propagates, the larger the loss 
of its memory. Accordingly, the spin filtration efficiency 
of the dsDNA is dominated by its end segment containing 
several base-pairs. 

To further verify the aforementioned point, we intro- 
duce point mutation in the dsDNA, where only one base- 
pair is modified and replaced by another. [36] Here, the 
point mutation is defined by switc hing the complemen- 
tary bases within single base-pair. [37] We focus on P s 
of the random dsDNA in Fig. Efa). The rd2 and rd3 
sequences are derived by introducing the point muta- 
tion in the rdl one.[37j One notes that P s is reduced 
more significantly if the mutation position is closer to 
the last base-pair of the rdl sequence. The largest P s 
is decreased from 42% for the rdl sequence to 25% and 
4.7% for the rd2 and rd3 sequences, respectively. (P s ) 
and (Gf ) are shown as a function of the mutation posi- 
tion L in Fig. E3d), where the average is obtained within 
the energy region [8.98,9.18]. It is clear that (P s ) does 
not change if the mutation occurs in the very beginning 
of the rdl sequence, and fluctuates more strongly if the 
mutation position becomes closer to its end. P s is very 
small if the point mutation takes place within the last 
three base-pairs, due to the identical sign between t\ n 
and t2n-3\ I n contrast, (Gf) fluctuates more obviously 
if the mutation occurs in the beginning of the sequence. 
And the fluctuation amplitude is more severe in the curve 
of (P s )-L than (G^)-L, indicating that the spin polariza- 
tion is much more sensitive to the modification of the 
base-pair in the dsDNA than the conductance. In this 
perspective, the spin transport along the dsDNA may be 
related to mutation detection in the biological processes 
and could be beneficial for DNA sequencing. [38[ 

Figure |4] shows the statistical properties of P s at fixed 
energy for the random and chr22-based dsDNA with 10 5 
samples. It clearly appears that P s can vary from 42% 
to negative, implying that the spin polarization direction 
of the charges transmitted through the dsDNA could be 
reversed by modifying its sequence. And one can see that 
many DNA molecules exhibit high P s . From a statistical 
perspective, the chr22-based dsDNA has more efficient 



spin filters than the random one. For instance, the num- 
ber of the dsDNA, of which P s is larger than 30% (20%), 
is 458 (1436) and 667 (2020) for the random dsDNA and 
the chr22-based one, respectively. This can be further 
demonstrated in the inset of Fig. 01 where one notices 
that the curve of the chr22-based dsDNA is always higher 
than that of the random one for P s > 1.3%. However, 
there are also many dsDNA with small P s at fixed en- 
ergy. This is attributed to the fact that: (1) P s depends 
on E that the energy region of large P s may differ from 
one sample to another [Figs.^a) and|3](b)]; (2) the elec- 
trons may be not polarized exactly along the helix axis 
for each dsDNA and the actual spin polarization could 
be larger. 

Finally, we study the spin polarization of the substi- 
tutional sequences of DNA molecules, of which the elec- 
tronic properties have been investigated previously, pjjl 
|40j Fig. [5] shows P s and G^ for several substitutional 
dsDNA. It is clear that both P s and G-j- are very large 
for these dsDNA. Therefore, besides homogeneous DNA 
molecules, other aperiodic DNA sequences can be also 
efficient spin filters with high spin polarization and con- 
ductance. 

In summary, we investigate the quantum spin trans- 
port through different dsDNA contacted by nonmagnetic 
leads. We find that the spin polarization strongly de- 
pends on the dsDNA sequence and is mainly determined 
by the end segment. Both natural and artificial dsDNA 
could be very efficient spin filters. Our results could mo- 
tivate further experimental studies on DNA spintronics. 
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